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I Abstract The propagation of nonhnear electron-acoustii 
waves (EAWs) in an unmagnetized cohisionless plasma 

I system consisting of a cold electron fluid, superthermal 
hot electrons and stationary ions is investigated. A 

[ reductive perturbation method is employed to obtain 
a modified Korteweg-de Vries (mKdV) equation for 
the first-order potential. The small amplitude electron- 
acoustic solitary wave, e.g., soliton and double layer 

I (DL) solutions are presented, and the effects of su- 
perthermal electrons on the nature of the solitons are 
also discussed. But the results shows that the weak 
stationary EA DLs cannot be supported by the present 
model. 

Keywords superthermal electrons; Electron- acoustic 
solitary waves; Double layers; mKdV equation 

1 Introduction 

The EAWs can either exist in a two temperature (cold 
and hot) electron plasma (Watanabe and Taniuti, 1977; 
Yu and Shukla, 1983) or in an electron-ion plasma 

[ with ions hotter than electrons (Fried and Gould, 
1961). EAWs are typically high-frequency (by compar- 
ison with the ion plasma frequency) , dispersive plasma 

> waves where the relatively cold inertial electrons oscil- 
late against a thermalized background of inertialess hot 
electrons which provide the necessary restoring force. 
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In the long-wavelength approximation, the dispersion 
relation of EAWs is given as a; = k (nco/nho)^^^ vth, 

1/2 

where vth = {nBTh/me) is the hot electron ther- 
mal speed, Kb is the Boltzmann constant, and n.co(^/io) 
are the cold (hot) equilibrium electrons densities. The 

1 /2 

phase speed Cg = {rico/nho) ^th of EAWs must be 
intermediate between cold and hot electron thermal ve- 
locities so that the Landau damping is avoided. In the 
Maxwellian plasmas, Gary and Tokar (1985) performed 
a parameter survey and found that the hot electron 
component constitutes a non-negligible fraction of the 
total electron density (more than ^20% ) for the exis- 
tence of the EAWs. And for the Lorentzian plasmas. 
Mace and Hellberg (1990) showed that the EAWs were 
usually strongly Landau damped by the hot electrons, 
unless rico/riho ^ 1. 

In the nonlinear wave studies, the propagation of 
solitary waves is important as it describes the char- 
acteristics of interaction between waves and plasmas. 
Solitary waves are localized nonlinear wave phenomena 
which arise due to a delicate balance between nonlinear- 
ity and dispersion. Among the best known paradigms 
used to investigate small-amplitude nonlinear wave be- 
havior are different versions of KdV equation (Washimi 
and Taniuti, 1966), or nonlinear Schrodinger equation 
(NLSE) (Hasegawa, 1975). Some form of reductive per- 
turbation technique is used to derive such equations. 
The KdV or mKdV describes the evolution of unmod- 
ulated wave, while the NLSE governs the dynamics of 
a modulated wave packet. In addition, for the arbi- 
trary amplitude solitary waves, the Sagdeev pseudopo- 
tential method (Sagdeev, 1966) is used too. Electro- 
static solitary structures are often observed in the space 
and laboratory plasma environment. EA soliton has 
been considered as one of the possible candidates for 
some of the observed solitary structures. Recently, the 
propagation of the linear as well as nonlinear EAWs 
has received a great deal of renewed interest not only 
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because the two electron temperature plasma is very 
common in laboratory experiments (Derfler and Simo- 
nen, 1969; Henry and Treguier, 1972) and in space 
(Dubouloz et al, 1991, 1993; Pottelette et al., 1999; 
Berthomicr ct al., 2000; Singh and Lakhina, 2001), but 
also because of the potential importance in interpreting 
electrostatic component of the broadband electrostatic 
noise (BEN) as being solitary EA structures observed 
in the cusp of the terrestrial magnetosphcre (Tokar and 
Gary, 1984; Singh and Lakhina, 2001), in the geomag- 
netic tail (Schriver and Ashour-Abdalla, 1989), in au- 
roral region (Dubouloz et al., 1991, 1993; Pottelette et 
al., 1999), in the numerical simulation (Lu, Wang and 
Don, 2005; Lu, Wang and Wang, 2005), and in labora- 
tory experiment (Lefebvre, et al., 2011), etc. 

On the other side, since Alfven and Carlqvist (1967) 
had suggested the current disruption theory for solar 
flare, the subject of DL (sometimes also called shock 
or kinks) has attracted great attention (Li, 1984; Liu, 
2010; reference therein). DLs occur naturally in a va- 
riety of space plasma environments. It turns out that 
DLs have the electrostatic potential and other relevant 
parameters monotonically changing from one value at 
one extreme to another at the other end, hence "kinks" . 
This is associated with adjacent positive and negative 
charge regions, which give rise to the name "double 
layers". Such DLs are more difficult to generate and 
require a fine tuning of the plasma parameters, hence 
a more complicated plasma compositions with enough 
leeway to obey the necessary constraints (Verheest, 
2006; Hellberg, 1992; Moslem, 2007). 

The study of nonlinear EAWs has been focused 
by many authors with different particle distribution, 
i.e.. Cairns distribution (Pakzad and Tribeche, 2010), 
Vortex-like distribution (Mamun and Shukla, 2002), 
q-nonextensive distribution (Gougam and Tribeche, 
2011), quantum plasmas (Masood and Mushtaq, 2008), 
et al. And it had been found that the particles distri- 
butions play a crucial role in characterizing the physics 
of nonlinear waves. 

Recently, the plasma with superthermal particles has 
gained much attention. Superthermal electrons are 
often observed in laboratory, space, and astrophysi- 
cal plasma environments, viz., the ionosphere, auroral 
zones, mesosphere, lower thermosphcre, etc (Picrrard, 
2010, and reference therein). The Kappa functions (Va- 
syliunas, 1968) characterized by the spectral index k 
are found to represent more suitably the particle's ve- 
locity distributions observed in number of space and as- 
trophysical environment. The Kappa function may re- 
cover the Maxwellian distribution in the limit of k — > oo 
and its mathematical characteristics and physical ori- 
gin have recently been addressed by Hau and Fu (2007, 



the references within). It is worth noting that some 
theoretical work focused on the effects of superthermal 
particles on different types of linear and nonlinear col- 
lective processes in plasmas. For instance, the linear 
properties of plasmas in the presence of a Kappa dis- 
tribution with excess superthermal particles have been 
investigated rather extensively (Summers and Thorne, 
1991; Mace and Hellberg, 1995). 

More recently, employing a Sagdecv pseudopotential 
method, the nonlinear arbitrary amplitude EAWs were 
studied, the weak stationary solitons and DLs were also 
given by expanding the Sagdeev potential in small am- 
plitude limit (Sahu, 2010; Younsi and Tribeche, 2010). 
In their model, the plasma systems were assumed con- 
sisting of cold fiuid electrons, superthermal hot elec- 
trons and stationary ions. The same procedure also 
was founded for the q-nonextensive distributed hot elec- 
trons plasma system (Gougam and Tribeche, 2011). 
Baboolal et al. (1991) have showed that, when ion- 
acoustic DLs in a plasma with negative ions are consid- 
ered, one must be especially careful to ensure that one's 
solutions meet the criteria for convergence of the origi- 
nal expansions. Verheest (1993) has arrived at a simi- 
lar conclusion considering DLs in dusty plasmas. Mace 
and Hellberg (1993) showed that EA DLs can not be 
supported in an infinite, homogeneous, unmagnctized 
and collisionlcss plasma system consisting of cool fluid 
ions, cold fluid electrons and hot Boltzmann distributed 
electrons based on the mKdV model. Some same con- 
clusions were also can be found (Hellberg, et al., 1992). 

Thus, with these ideas in mind, we re-investigate 
here the small but flnite amplitude solitary structures in 
such two-electron-tcmperature plasma with superther- 
mal electron based on the mKdV model. One of our 
objective here is to study the nonlinear effects of super- 
thermal distribution of hot electrons on the nature of 
the small amplitude solitary waves. Another one is to 
show whether the DL solution exist or not. This pa- 
per is organized as follows: in Section 2, the basic set 
of equations is introduced. In Section 3, we derive the 
mKdV equation, and the solutions of both solitons and 
DLs are given. Finally, some conclusions and discus- 
sions are given in Section 4. 

2 Basic equation 

We consider a homogeneous system of an unmagnc- 
tized collisionless plasma consisting of a cold electron 
fluid, and superthermal hot electrons obeying a Kappa 
distribution, and ions. It is well known to us that 
the linear spectrum of EAWs extends only up to the 
plasma frequency of the cold electron population, ojpc = 
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y 47rncoe^/me, so at this high frequency, the ion popu- 
lation plays the role of a neutralizing background. For 
a small — but finite amplitude waves propagate one di- 
mensionally, such system is governed by the following 
normalized equations: 

Ql^C + ^ M = 0, 
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and the super-thermal hot electrons density is given 
by (Saini, 2009; Rios, 2010): 
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(2) 



where ndnh) are cold (hot) electrons normalized densi- 
ties to the total unperturbed density ng = rico + nho, Uc 
is the velocity of cold electrons normalized to the EA 
velocity Cg, (p is the electrostatic potential to KsTh/e, 
respectively. Time and space variables are normalized, 
respectively, to the inverse of cold electron plasma fre- 
quency uip^ and the hot electron Debye length Ad/i = 

^ KsTfi/ 4:Trnhoe'^ ■ The parameter a = nho/uco- It 
should be noted that the superthermal hot electron 
component constitutes a non-negligible fraction of the 
total electron density, as mentioned in the introduction. 



3 The mKdV equation and the solution 

Nonlinear electron-acoustic waves are governed by the 
full set of Eqs.((I]) and To derive the mKdV equa- 
tion describing the behavior of the system for longer 
times and small but finite amplitude EAWs, we employ 
the familiar reductive perturbation technique (Washimi 
and Taniuti, 1966; Mace and Hellberg, 1993). From the 
usual considerations of the small-wavenumber disper- 
sion relation for electron-acoustic waves, we introduce 
the slow stretched coordinates: ^ = e {x — Xt) , t = e^t, 
where £ is a small dimensionless expansion parameter 
and A is the wave speed normalized by Cg. All physical 
quantities appearing in Eq. ([T]) are expanded as a power 
series in e about their equilibrium values as: 
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We impose the boundary conditions that as |^| — ?> 
oo, ric — l,Wc = 0, = 0. Substituting Eq.Q into the 
system of Eq.(IT]), and equating coefficients of different 
powers of e, from the lowest-order equations in e, the 
following results are obtained: 
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Eqs.® and Q as well as Poisson equation give the 
linear dispersion relation as 
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From the next order of e, with the aid of equation 
we obtain the following equations: 
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Poisson equation with the help of Eq. ([2]) at O [e^ 
yields 



_ ^,3^ (2k-1)(2. + 1) 

2A4 2 (2k -3)' 

Since we assume that (p^^^ ^ 0, it follows that at least 
\A\ ^ 0(e). Hence it should be included in the next 
higher order, i.e., O (e^) of Poisson equation. Collect- 
ing the third order of e, that in O (e'^) , from equations 
(dl) and ([2]) using we obtain 
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Eliminate the third-order perturbed quantities n, 
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u'c \ and 'p^'^\ we obtain the following mKdV equation 
for the first-order perturbed potential: 

^^,,3a4.^ + «3|.' + |,.^0, ,11) 

with the coefficients read as 

(2k- 1) (2k 1)' 
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111 Ea. ([TT]) . (y9 is used in place of (^'^^ for brevity. Let 
us introduce the variable, ?? = $ — Mqt, where r\ is 
the transformed coordinates with respect to a frame 
moving with velocity Mq, for the steady-state solution 
of the mKdV equation pT|) . By using the boundary 
conditions iy9 and dLp/drj — > at \rj\ — > oo, we 
obtain 



dip 
drj 



+ y(^) = o, 



(14) 



where V [ip) is the Sagdeev pseudopotential (Sagdeev, 
1966), reads as 



V ((/?) ^ Ai^p^ + A2(p^ + Asif'^ 



with Ai = — A/o/A"^. It should be noted that when 
Mo < A, the resuh of Sahu (2010), i.e, in which Eq.(18) 
is recovered for the small amplitude solitary waves by 
Sagdeev pseudopotential approach. 

3.1 Solitons solution 

If one neglects the term corresponding to ip'^ in the 
Sagdeev pseudopotential p^ . the solution of Ea. ([T4| 
is 
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where the soliton peak amplitude pm and width A 
are given by pm = —AxjAi, A = ij \J~^2A[. Solu- 
tion represents a small-amplitude stationary EAWs 
provided Ai < or Mq > 0, which means that the soli- 
tary waves are supersonic, in agreement with the large 
amplitude case. It is clear that the nature of the soli- 
tary waves, i.e., whether the system will support com- 
pressive or rarefactive solitary waves, depends on the 
sign of A2. If A2 is positive (negative) a compressive 
(rarefactive) solitary wave exists. In our present case, 
< since k > 3/2 for a physically realistic thermal 
speed (Summers and Thorne, 1991; Mace and Hellberg, 
1995), so that here we would have a rarefactive soliton. 

3.2 Double layer solution 

For the DL solution, the Sagdeev potential must satisfy 
the conditions (Mace and Hellberg, 1993), (a): V {ip) = 
at (/J = and ^) = (b): V {p) = at = 
and = '■Pm ; and (c): V" {ip) < at = and 



p) = ipm- Applying the boundary conditions (a) and 
(b) into Eq. (fTSj) . we obtain 
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2A3' " 4^3 
then Sagdeev pseudopotential (|T5|) can be written as 

V{p>)= AiP>''[^ra-^f ■ (18) 

The DL solution is then given by (if it exists) 
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with A = -^Z— 8/A3 1 \p)m\ represents the width of the 
DL provided A^ < 0, i.e.. 



(15) ^ - 
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For a DL exists, the condition ([20| must be fulfilled, 
then we have that the hot e lectron concentration a must 



satisfy a < 



(2k+1)(2k+3) 



It is indicates that a < 1 in 



15(2k-1)^ 

the both small and large value of superthermal index 
K. At such large cold electron density ratio (n-co/'^ho > 
1), linear Landau damping by the cold electrons would 
become appreciable (Gary and Tokar, 1985; Mace and 
Hellberg, 1990). Furthermore, we note that for this case 
the mKdV model leads to a value of |0,„| ^ 1 in the 
small superthermal index case and > 1 in the large 
superthermal index case, which is shown in the Fig.l. 
It clearly does not satisfy the expanding requirement 
\4>m\ £• From the aforementioned considerations, we 
are let to conclude that Eq .(fT9| is not a valid solution of 
the mKdV equation for electron-acoustic waves in the 
present model. So we can draw the conclusion that the 
present model may not sustain a double layer structure. 
Such a point has also been studied extensively by Mace 
and Hellberg (1993) for the two-temperature electrons 
Maxwellian plasmas. Indeed, a systematic investigation 
should be be undertook to verify this conclusion, but it 
beyond the scope of present paper. 



4 Remarks and Conclusion 

To investigate the existence regions and nature of the 
EA solitary structures, we have done numerical calcu- 
lations for different set of parameters. 

Figs. 2 and 3 present the effects of the superthermal 
index k and hot electron concentration (a) on the am- 
plitude and width of the solitons in the slightly super- 
sonic point. It can be seen from Figs. 2 and 3 that the 
decreasing the superthermal index k and increasing in 
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hot electron concentration (a) will decrease the ampli- 
tude of soliton. 

In summary, we have studied the nonlinear EAWs in 
an unmagnetized, collisionlcss plasma consisting of cold 
electrons, superthermal hot electrons, and stationary 
ions. A reductive perturbation method has been used 
to get the mKdV equation which describes the dynam- 
ics of solitons and DLs. The effects of super-thermal 
index k and concentration a of the hot electrons to the 
cold electrons on the nature of the solitons are also dis- 
cussed. It is seen that the effects of the superthermal 
hot electrons have a very significant role on the ampli- 
tude and width of the weak amplitude EA solitons. It is 
also found that the small amplitude DLs cease to exist 
in our present model. Considering the wide relevance of 
nonlinear oscillations, we stress that the results of the 
present investigation should be useful in understanding 
the nonlinear features of localized electrostatic acoustic 
structures in different regions of the laboratory experi- 
ments and space environments. 




Fig. 2 Small amplitude soliton for different 3 for 

solid line, k = 6 for dotted line, and k = 15 for dashed line, 
where the other parameters are a — 0.5 and Mo = 0.1 A. 



Acknowledgements We thank Professor X Q Li 
for the helpful suggestions and encouragement. This 
work is supported by the National Natural Science 
Foundation of China under the grant Nos. 10963002 
and International S&T Cooperation Program of China 
(2009DFA02320) and the National Basic Research Pro- 
gram of China (973 Program)(No. 2010CB635112). 

The authors are grateful to the referee for his com- 
ments for the improvement of this paper. 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



Fig. 1 Plot of EA DLs amplitude Ifml vs k, where a 
0.3, 0.25, 0.2 for solid, dashed, dotted line, respectively 




Fig. 3 Small amplitude soliton for different a, i.e, a — 0.3 
for solid line, a = 0.5 for dotted line, and a = 1.0 for dashed 
line, where the other parameters are n — 3 and Mq = O.IA. 
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